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By combining the recent WMAP measurements of the cosmic microwave background anisotropies 
and the results of the recent luminosity distance measurements to type-la supernovae, we find that 
the normalization of the matter power spectrum on cluster scales, erg, can be used to discriminate 
between dynamical models of dark energy (quintessence models) and a conventional cosmological 
constant model (ACDM). 
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I. INTRODUCTION 
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The WMAP satellite measurements of the Cosmic Mi- 
crowave Background anisotropies P| have provided accu- 
rate determinations of many of the fundamental cosmo- 
logical parameters. When combined with other data sets 
such as the luminosity distance to type-la supernovae or 
large scale structure (LSS) data they rein- 

force the need for an exotic form of dark energy, which 
is characterized by a negative pressure and is responsible 
for the observed accelerated expansion of the universe. 
There are two main scenarios used to explain the na- 
ture of the dark energy, a time independent cosmological 
constant A and quintessence, which involves an evolving 
scalar field Q U [S, 9] . Previous tests of quintessence 
with pre- WMAP CMB data [ll|ll|l2, have led to con- 
straints on the value of the dark energy equation of state 
parameter, wq ^ —0.7 with the cosmological constant 
value, WA — —1 being the best fit. Nevertheless a dy- 
namical form of dark energy is not excluded. Specifi- 
cally the detection of a time variation in this parameter 
would be of immense importance as it would rule out 
a simple cosmological constant scenario. When parame- 
terising quintessence models we do not want to assume 
simply a constant equation of state wq since this intro- 
duces a systematic bias in the analysis of cosmological 
distance measurements jl3i |. with the effect of favouring 
larger negative values of wq if the dark energy is time 
dependent. For instance it is possible that claims for a 
'phantom' component, where wq < —1 0, Q| are en- 
tirely caused by this effect. Moreover assuming wq con- 
stant underestimates the contribution of the dark energy 
perturbations (which are a specific feature of scalar field 
models) on the evolution of the gravitational potentials 
and consequently the effect on the CMB power spectrum 
|15|. In this paper we deliberately do not assume wq to 
be constant, rather we focus on the relation between a 
dynamical dark energy component and the normalisation 
of the dark matter power spectrum on cluster scales, erg- 
We also discuss the age of the universe, to, and show how 
the new data sets undermine its use for distinguishing be- 



ll. METHOD AND DATA 

In this analysis, rather than considering a specific 
scalar field model, we allow for a time dependence of 
the equation of state parameter w q. Several formulae 
have been proposed in the literature [19. Il7j all with lim- 
ited applicability. In ll^l a form for wq (z) was suggested 
which is valid at all redshifts and parameterises the equa- 
tion of state in terms of five parameters, which specify 
the value of the equation of state parameter today Wq, 
and during the matter/radiation eras Wq/wq] the scale 
factor a™ where the equation of state changes from Wq 
to Wq and the width of the transition A. Since Big-Bang 
Nucleosynthesis bounds limit the amount of dark energy 
to be negligible during the radiation dominated era, with- 
out loss of generality we can further reduce our parameter 
space by setting Wq — Wq in Eq. (4) of Ref . . The 
parameters given by the vector Wq — (wg, Wq , a™. A) 
can account for most of the dark energy models proposed 
in the literature. For instance quintessence models char- 
acterized by a slowly varying equation of state, such as 
supergravity inspired models [l9| , correspond to a region 
of our parameter space for which < a™/ A < 1, while 
rapidly varying models, such as the two exponential po- 
tential case 20], correspond to a™ /A > 1. Models with a 
simple constant equation of state are given by Wq ~ Wq . 
The cosmological constant case is also included and cor- 
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responds to the following cases: 

lyO ^ < Q.f with a™/ A > 1. Assum- 

ing a flat geometry we perform a likelihood analysis of 
the WMAP data to constrain dark energy models speci- 
fied by the vector Wq and the cosmological parameters 
Wc = {flQ,i^bh^ ,h,ns,T,As) which are the dark en- 
ergy density, the baryon density, the Hubble parameter, 
the scalar spectral index, the optical depth and the over- 
all amplitude of the scalar fluctuations respectively. We 
have modified a version of the CMBFAST code to 
include the dark energy perturbation equations in terms 
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of the time derivatives of the equation of state [23|- In 
order to break the geometric degeneracy between Wg , fig 
and h, we use the most recent compilation of supernova 
data of 4] in addition to the WMAP TT and TE power 
spectrum data. We evaluate the likelihood of CMB data 
with the help of the software provided by the WMAP 
team |2,Sj . The important point which we want to stress 
is that we are able to treat both data sets, (WMAP and 
SN-Ia) without making any prior assumptions concern- 
ing the underlying cosmological model, in order to be as 
conservative as possible and to evade potential problems 
with issues like relative normalisations and bias. We re- 
strict our analysis to dark energy models that satisfy the 
null dominant energy condition and Wq , Wq > — 1 and 
following the analysis by the WMAP team, we use the 
prior T < 0.3 in order to prevent Qb from taking unphys- 
ically high values. 



III. RESULTS 

The WMAP CMB data constrains the cosmological 
parameters Wc in a range of values consistent with the 
results of previous analysis such as 0, |3j Ull • In partic- 
ular we find the scalar spectral index ns = 1.00 ± 0.04, 
the physical baryon density ilfc/i^ = 0.0234 ± 0.0014 and 
the optical depth r = 0.17 ± 0.06. As mentioned above, 
in order to break the degeneracy between Wq, flq and 
h, we combine the CMB data with the SN-Ia luminos- 
ity distance measurements. This allows us to constrain 
the Hubble constant to be /i = 0.68 ± 0.03, in agree- 
ment with the HST value [2^, the dark energy density 
fig = 0.72 ± 0.04 (all limits so far at la) and the present 
value of the equation of state Wq < —0.82 (at 95% CL). 
It is important to stress that the the addition of the dark 
energy parameters Wq does not introduce any new de- 
generacies with the other parameters. This is clear from 
the fact that the constraints on Wc are in agreement 
with other previous data analyses. Figure ^ shows the 
marginalised one-dimensional likelihoods for ACDM and 
the dynamic dark energy models. We will defer a detailed 
discussion of these results to a later paper, and in this 
paper concentrate on the use of dark matter clustering 
as a probe of quintessence models. 

In general we expect dark energy to affect the value 
of (Tg because it can lead to a different expansion history 
of the universe 12^ . However, in it was shown that 
different dark energy models leave particular imprints on 
the large angular scales of the CMB anisotropy power 
spectrum through the integrated Sachs- Wolfe (ISW) ef- 
fect. The excess of power produced by the ISW at low 
multipoles affects the normalization of the matter power 
spectrum 28]. For instance models with a fast late time 
transition in the equation of state produce a larger ISW 
effect than a pure cosmological constant scenario. As a 
consequence they require a smaller amplitude of primor- 
dial fluctuations in order to match the observed CMB 
spectrum. In this case the predicted value of erg will be 




0.2 0.3 0.4 




60 65 70 75 80 




J " " 1 

~ / 

: / 

" ' / 
1 1 


" " 1 ' " ' I 

X\ 

' \ \ 
\ \ 


~ ' 1 
// 

~ // 

'■_ il 


\ \ - 


\ y 


\ \| 


~ '/ 

i/r . . . 1 


\ - 
— 1 . . .N 



0.1 0.2 0.3 



FIG. 1: Marginalized likelihoods for the various cosmological 
parameters in the ACDM scenario (red dashed curve) and 
including the QCDM models (black curve). The results agree 
very well with each other. 



smaller than in the ACDM model. This specific class of 
models has already been investigated using pre- WMAP 
data [l^ |2^ . but the results underestimated the opti- 
cal depth subsequently found by WMAP, leading to an 
over-estimation of the power on small angular scales. It is 
only with the release of the first year of WMAP data that 
through one CMB data set, we can link the anisotropics 
on large and small angular scales. This is an exciting 
feature of the data, as it allows us to properly assess the 
effects of ISW and the normalization of the matter power 
spectrum. In figure |21 we plot the two dimensional likeli- 
hood contours in the f2„i — cg plane. The filled contours 
correspond to 1 and 2ct values for the dark energy models 
spanned by Wq, while the solid curves correspond to the 
ACDM case. As expected from the above discussion, we 
note that A models have systematically higher values of 
CTg than models with a time varying equation of state. 

It seems clear that a CMB independent estimate of the 
value of (Tg would be able to distinguish between a ACDM 
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FIG. 2: Marginalized 68% and 95% confidence contours for 
quintessence (filled contours) and ACDM models (solid lines). 
ACDM has a systematically higher value of erg, and a slightly 
higher value of $7^. 



and dynamical equation of state model. For instance val- 
ues of (78 < 0.7 would be rejected at 2cr in the ACDM 
case. More specifically we find that the value of a% can 
discriminate between different dark energy models. This 
can be seen in figures O ^ and |S1 which are the main re- 
sult of this paper. In Fig. O we plot the average value 
of (Ts as a function of a™ and Wg, where the average is 
taken over all models in our chain which exhibit a rapid 
transition (defined here as w'q > —0.2 and A < 0.1). 
A ACDM model corresponds to a™ and Wq = — 1. 
The average value of ag in this point is 0.9. As we move 
away from the ACDM corner, the average erg decreases 
monotonically, as seen by the contours. To assess the 
usefulness of erg for distinguishing between models given 
todays data, we also plot two 68% confidence regions, 
one for models with cg > 0.9 (lighter gray) and one with 
(Tg < 0-6 (darker gray). Clearly, if we restrict ourselves to 
models with a high value of ug, we favour a ACDM-like 
behaviour of the dark energy. In the opposite case, we 
find a™ ~ 0.3. Together with the fast-transition condi- 
tions given above, this means that these models have an 
equation of state w{z > 2) :$> —1, and we would exclude 
the case p = —p at over 95% CL. As we marginalise over 
all other parameters, we see that no degeneracies spoil 
this result. 

As a complementary view, we can plot a™ and Wq for 
fast-transition models (without the condition on Wq), see 
Fig.m The data requires that < -0.8 and so ACDM 
models occupy the region defined by either a™ —>■ (in 
which case the equation of state is independent of Wq) 
or Wq — > — 1 (and thus w{z) « — 1 without transition), 
which again coincides with the high-ug models. Models 
with (Tg < 0.6 on the other hand require both a™ '~ 0.3 
and w7i ~ -0.7 at 68% CL. 



Fig-Eis the corresponding figure for dark energy mod- 
els with a slowly varying equation of state (0 < a™/ A < 
0.8). In this case the relevant parameters are Wq and 
Wq, and the ACDM models are now at Wq = Wq = — 1. 
Again, erg decreases rapidly as we move away from that 
corner. We show once more the Icr regions for models 
with CTg > 0.9 (lighter gray) and with trg < 0.6 (darker 
gray). Models with a high value of erg are again clus- 
tered around the ACDM region, and those with a low 
clustering amplitude require w ^ —1 at high redshift. 

We expect these regions to shrink as the cosmologi- 
cal parameters become more constrained by future data, 
which will improve the impact of clustering as a probe 
of the time dependence of the dark energy. This is our 
main result, and it means that, given a precise measure- 
ment of CTg, we can impose strong limits not only on 
the value of w today, but also at earlier times. Even 
if Wq « — 1 today, we are able to probe its behaviour 
at higher redshift and to either exclude ACDM or sig- 
nificantly constrain quintessence type models. Although 
especially slowly varying models cannot be ruled out as 
they can approximate the behaviour of a true cosmologi- 
cal constant arbitrarily closely, these models become less 
and less attractive as they start to require the same fine 
tuning as A itself. 
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FIG. 3: The average as as a function of Wq and a™ for 
models with a rapid transition in wq (numbered lines). We 
also show the 68% confidence regions for models with erg < 0.6 
(dark gray) and erg > 0.9 (light gray). 

Why are we using ug as a variable as opposed to simply 
choosing one of the many published measured values of 
(Tg? First, the published data shows a large spread of val- 
ues 's^l, so that our conclusions would strongly depend 
on the choice of data sets. Secondly, the measurements 
also depend in general on the dark energy parameters and 
the results quoted are only valid for ACDM models. For 
example, this is the case for the large scale structure re- 
sults, which implicitly assume a ACDM model when pass- 
ing from redshift space to real space, and for weak lens- 
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FIG. 4: The average as a function of lUg and a™ for 
models with a rapid transition in wq (numbered lines). We 
also show the 68% confidence regions for models with erg < 0.6 
(dark gray) and as > 0.9 (light gray). 
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FIG. 5: The average erg as a function of uiq and for 
models with a smoothly varying wq (numbered lines). We 
also show the 68% confidence regions for models with as < 0.6 
(dark gray) and ag > 0.9 (light gray). 



ing measurements. In the second case, the dependence 
on the dark energy characteristics is strong enough that 
it can be used to constrain the evolution of the equation 
of state [sj. As an illustration, we can assume that the 
clustering results deduced from velocity fields in Ref. 1^3 
are unaffected by the details of the dark energy evolution. 
As a rough approximation to their PSCz results, we set 
(78 « (1.13 ± 0.05)(fl„/0.3)° ^ In this case, the con- 
straints on quintessence models become much stronger, 
e.g. < -0.9 at 95% CL. On the other hand, if fu- 
ture precision measurements converge on cts ^5 0.7 then 
ACDM is ruled out at high significance. 

Moreover, erg is linked to the amplitude of the mat- 



ter power spectrum P{k) on small scales. To measure a 
possible running of the scalar spectral index, dng/dlogfc, 
in inflationary models, it is necessary to combine CMB 
data on large scales with P{k) on small scales. Since 
quintessence models can change the amount of cluster- 
ing on small scales with respect to a ACDM model, it is 
possible for them to mimic the effect of such a running. 
This possibility should be kept in mind when constrain- 
in^ models through the combination of different data sets 




FIG. 6; Marginalized 68% and 95% confidence contours for 
quintessence (filled contours) and ACDM models (solid lines). 

Another observable which has been studied in this con- 
text is the age of the universe, to H^, 13 1 which is in 
general also a function of the dark energy parameters 
Wq . An independent measurement of (for which the 
WMAP limit does not qualify, as it explicitly assumes 
ACDM) can thus be used to set limits on the equation 
of state. Since the luminosity distance and to pos- 
sess a similar dependence on the Hubble rate, the SN-Ia 
data, which probe about two-thirds of the age of universe, 
can provide tight constraints on tneven for generic dark 
energy models. For example in [J| considering ACDM 
cosmologies, the authors obtain Hoto — 0.96 ± 0.04. The 
limit is also valid for quintessence, as we find Hoto = 
0.96 ± 0.03 for the combination of CMB and SN-Ia data. 
This constraint, together with the remaining slight de- 
generacy in Hq which leads to lower values of the Hub- 
ble constant as we move away from the ACDM models, 
means that the allowed quintessence models are older 
than the those with a cosmological constant, as we can 
see in Fig. |B| The marginalised age of quintessence uni- 
verses is to — 13.8 ± 0.3 Gyr, while in the ACDM case 
to — 13.55 ± 0.26. Clearly, it will be difficult to use to 
to disentangle different models until the uncertainty in 
the cosmological parameters is further reduced. But if 
we were to find a lower limit on the age of the universe 
which is too high for ACDM, we could potentially inter- 
pret it to be a sign of quintessence. 
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IV. CONCLUSIONS 

In this paper we have demonstrated how, by combin- 
ing WMAP and SN-Ia data, it is possible to use the nor- 
mahsation of the dark energy power spectrum on cluster 
scales, (Tg, to discriminate between dynamical models of 
dark energy (quintessence models) and a conventional 
cosmological constant model (ACDM). In particular we 
have shown for the first time that a CMB independent 
measurement of erg allows us to constrain the parame- 
ters describing the evolution of the dark energy equation 
of state. For instance, we found that standard ACDM 
is ruled out at over 95% CL (compared to a time depen- 
dent dark energy component) if erg < 0.7. This constraint 
can be relaxed by going beyond the standard model, i.e. 
introducing very massive neutrinos or a running of the 
spectral index [S^. However, we expect improved data 



to lead to stronger limits in the near future. We have also 
briefly discussed the use of the age of the universe to as a 
way of constraining dark energy models, and shown that 
by itself it does not discriminate between quintessence 
and ACDM models, although coupled with ag, it may 
act as a useful cross check. 
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